Exactly solvable model of three interacting particles in an external magnetic field 
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The quantum mechanical problem of three identical particles, moving in a plane and interacting 
pairwise via a spring potential, is solved exactly in the presence of a magnetic field. Calculations 
of the pair-correlation function, mean distance and the cluster area show a quantization of these 
parameters. Especially the pair-correlation function exhibits a certain number of maxima given by 
a quantum number. We obtain Jastrow pre-factors which lead to an exchange correlation hole of 
liquid type, even in the presence of the attractive interaction between the identical electrons. 



I. INTRODUCTION 

The Laughlin wave function [1,2] was introduced to de- 
scribe the ground state of a two-dimensional (2D) long- 
range correlated electron gas in the fractional quantum 
Hall regime. The strong short-distance repulsion between 
any two particles is described by the Jastrow prefactor in 
the wave function. The same Jastrow prefactor appears 
also in the exact wave functions for N particles interact- 
ing pairwise through a potential of the form l/(xi — Xj) 2 , 
first obtained in one dimension by Calogero [3,4], and 
later extended to two and higher dimensions [5-7] . These 
wave functions lead to the Wigner-Dyson distribution 
function [8] for N particles. 

A natural question therefore is whether the origin of 
this Jastrow prefactor is due to special features of the 
interactions mentioned above, or it is a generic property 
of models with short-range interactions. 

The recent advances in the theory of disordered sys- 
tems show that the distribution function of physical pa- 
rameters of the mesoscopic system, such as the density of 
states (DOS) or the conductance, is represented by the 
Wigner-Dyson function in the weak localization regime 
(see, e.g., [9] for a review). Physically, this describes a 
repulsion between the energy levels. Averaging the par- 
tition function of a disordered system over the Gaussian- 
distributed one-particle impurity potentials results in an 
effective Hamiltonian with short-range interactions. This 
clearly shows that level repulsion is not a property of just 
few special interaction potentials. 

In a recent paper [10], a new mechanism for the for- 
mation of three-particle clusters in Si-MOSFET struc- 
tures and GaAs/AlGaAs heterojunctions has been stud- 
ied. The exchange type interactions between 2D band 
electrons in the inversion layer and charged impurities 
in the oxide of the MOSFET can lead to an effective 
three-particle attractive interaction. A qualitative anal- 
ysis shows that the ground state of the three particles be- 



comes unstable with respect to arbitrary small attractive 
interactions between them. An attraction between three 
particles leads to the formation of bound states with neg- 
ative energies [10]. Although a weak attractive three- 
particle interaction does not produce a transition to a 
liquid state of the electron gas as described in Laughlin's 
approach, the ground state with three-particle clustering 
could be energetically favored in the fractional quantum 
Hall regime. This is a consequence of the occurrence of 
an effective pairing [10], which further reduces the ground 
state energy. 

It is interesting to study the internal structure of a 
particular cluster with three particles, in order to find 
out if its wave function contains a Jastrow prefactor. This 
would prevent the particles from approaching each other 
even in the attractive case. In this paper we study this 
question by modeling a generic three-particle attractive 
potential by a pairwise spring-like interaction. We show 
that this choice of potential allows for an exact solution 
even in the presence of a constant external magnetic field. 

Exact results obtained for solvable models are im- 
portant to understand general features of three-particle 
problems and to test approximate solutions, such as the 
ones obtained from the widely used Fadeev's equations 
[11]. Very few exact solutions of three-particle problems 
have been found, mostly for one-dimensional [12-14] or 
for three-dimensional [15-17] spinless particles. In this 
paper we present an exact solution for a model of parti- 
cles with spin. 

The topology of the three-particle cluster may be either 
of a string form, when two particles with opposite spin 
are placed in the same spatial point, or of a triangle form 
when all three particles lie far from each others. In the 
former case the number of the spin configurations is equal 
to two (doublet states) with total spin S = 1/2. Apart 
from this doublet states, a triangular clustering of quartet 
states may be also realized. This would correspond to 
four different symmetric spin wave functions with a total 
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spin 5 = 3/2. In a strong magnetic field this quartet 
states are the only relevant ones. 



II. THE MODEL 



*(ri,r 2 ,r 3 ) = ip(r)ip(^)tp(r)). 



(7) 



Schrodinger's equation Hip — Eip with the Hamiltonian 
given by one of the Eqs. (4)- (6) is easily solved in a polar 
coordinate system where, e.g., Eq. (6) is written as 



The Hamiltonian of the model is 

~ h 2 A / d . e 

H = > i— A a 

2m* ^ \dr a he 

a=l 

+ |[(r 1 -r 2 ) 2 + (r 1 -r 3 ) 2 + (r 2 -r 3 ) 2 ], (1) 

where r a (a = 1,2,3) are the two-dimensional (2D) po- 
sition vectors of electrons with effective mass m* . The 
spring strength is chosen as n — m*oj 2 /3 where ui is 
the specific frequency of the electron's relative vibration. 
A a = -§{— y,x,0} is the symmetric-gauge vector poten- 
tial at the position of the a-th particle for a magnetic field 
B = {0, 0, B/2} perpendicular to the electronic inversion 
layer. 

In order to diagonalize the Hamiltonian we introduce 
the new coordinates {r, 

1 , 

? = 7! (ri - r2) ' 



ri +r 2 



1-3 



(2) 



where r is the center-of-mass coordinate, £ is the relative 
coordinate of the particles 1 and 2 and rj is the relative 
coordinate of the third particle with respect to the center- 
of-mass of the particles 1 and 2. In the new coordinates, 
the Hamiltonian is decoupled, 



H — H r + + H v , 



(3) 
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(V x V„) z . 
(6) 



The last terms in the Hamiltonians given by Eqs. (4)-(6) 
are proportional to the angular momentum and 1% = j§ 
is the magnetic length. 

The wave function W of the initial Hamiltonian, Eq. 
(1), can be presented as the product of three wave func- 
tions tp which are eigenf unctions of Eqs. (4)- (6) 



H P = - 



h 2 / d 2 id_ l o> 2 \ . h 2 d 

2m*\dp 2 pdp p 2 dip 2 J 2m*l B dtp 



+ 
3k n 
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(8) 



and p and <p are the radial and angular variables in the 
polar coordinate systems. 

The Hamiltonian (8) commutes with the angular mo- 
mentum operator L z = —ifr-§^, [H p ,L z ] = 0. Therefore 

the wave function ij>(p, ip) of the Hamiltonian H p can be 
chosen to be of the form 



1> m {p,<p)=e im *R(p), 



(9) 



where m is the angular momentum for a 2D system and 
to = ±1,±2,±3.... 

The asymptotic behavior of the radial wave function 
R(p) can be separated as 



R(z) = z S ? L e-t'f(z), 



(10) 



where z = \ 1 1 + -t 2 - is a dimensionless coordinate, 

ii b y u b 

and lob = -Nr- is the cyclotron frequency. 

The function f(z) satisfies the confluent hypergeomet- 
ric equation 



f(z) = F{-( 



V 71 + 4"o/"b 



\{\m 



+ l)),l + |mU} 
(11) 



where (3 = — tt ^ the dimensionless energy spectrum. 
Convergence requires the first coefficient of the confluent 
hypergeometric function to be a negative integer. The 
function then reduces to the generalized Laguerre poly- 
nomials and gives the following energy spectrum 



E ni ,mi = fl^B 




\mj\ + l hoj B 

o ) + —^- m h (12) 



where i = 1,2,3 denotes the number of non- interacting 
quasi-particles introduced by means of a normal coor- 
dinate transformation. The total energy of a three- 
particle cluster is expressed as a sum of the energies 
of three non-interacting quasi-particles. Redefining the 
initial quantum numbers as no = n\ + ^(\mi \ + mi), 
n = 2(n 2 + n 3 ) + |m 2 | + |m 3 | and m = m 2 + m 3 , the 
total energy becomes dependent on three exact quantum 
numbers, 
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E(no, n, m) = huu B {no + -) 



+ - 



'l + ^(„ + 2) (13) 



The case corresponding to three free particles in an ex- 
ternal magnetic field is given by the two exact quantum 
numbers n = J2i=i( n i + l m i|/2) an d m — 5Z i=1 m i- The 
normalized wave function of the i-th quasi-particle can 
be written in the form 
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(14) 



In order to construct the total wave function for a three 
particles cluster, the spatial wave functions correspond- 
ing to Eqs. (4)-(6) must be multiplied by the appropriate 
spin wave functions and the final expression must be an- 
tisymmetric. 



III. WAVE FUNCTIONS AND CORRELATORS 



All spin wave functions of the quartet states \S = 



|,s z ), corresponding to the total spin number S 
and its z-component s z , are symmetric 

|3/2, 3/2) - a ia2 a 3 

|3/2, 1/2) = -^={ ai a 2 (3 3 + ai (3 2 a 3 + (3 ia2 a 3 } 

|3/2, -1/2) = -^={ ai p 2 p 3 + /3i/3 2 a 3 + p ia2 p 3 } 
|3/2, -3/2)= 



= 3/2 

(15) 
(16) 

(17) 
(18) 



where a* and 0i are the spinors of the i-th electron cor- 
responding to spin up and down, respectively. Therefore, 
the antisymmetric orbital part of the wave function reads 

* (ri,ra,r s ) = C* ni>mi (r; w = 0) 

X{*n 2 ,m 2 (0*ns,ms(»7) - *n 2 ,m 2 (-0*n 3 ,m, (v) 
+*n 2 ,m 2 (0* (»?) - *n 2 ,m 2 (-0*n3,m3(»7) 

113,113 ft) "^n^ro^ 0^"3>"i3 

where C is a normalization constant and we have used 
the notation 
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FIG. 1. The normalized pair correlation function <?(£*) ver- 
sus the dimensionless distance between two particles for 
712 = 1 and different values of the quantum numbers ri3 , rri3 . 
The dotted lines are the pair correlation function only from 
the diagonal elements of |*| 2 according to (25). 

Notice that the inversion of the coordinates in polar 
basis leads to p — > p and <p — > <p + tt. Under this transfor- 
mation, the wave function acquires a factor (— l) m , which 
defines its parity. It is easy to see that ^(ri, r 2 , r 3 ) = 
for m 2 = even integer. The privileged dependence of the 
wave function on m 2 is originated by the asymmetry in 
the new coordinates (2). For m 2 = odd integer, we get 

* Q (ri,r 2 ,r 3 ) = Ctf nilTni (r; w c = 0) 

(0*™ 3 ,m 3 (f?) 
+ *n 2 ,m 2 (0* (v) 
+ *n 2 ,m 2 (0* n 3 ,m 3 

(21) 

In high magnetic fields the polarized spin structure is 
characterized by this quartet states since all states are 
aligned. The valuable feature of the spin-polarized quar- 
tet states is that they permit only odd values of m 2 
and exclude even ones, which is the situation observed 
in the fractional quantum Hall experiments. The wave 
function ^ given by Eq. (21) has the property that 
^^(ri, r 2 , r 3 )| ri ^ rj = for arbitrary chosen i ^ j = 
1, 2, 3 . This corresponds to the behavior of the Jastrow 
prefactor /(r^ — rj) in Laughlin's ground state wave func- 
tion [1,2]. 
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since m 2 is odd. The regime of £* < 1 can be realized, 
according to Eq.(23), for an arbitrary finite distance be- 
tween particles by sufficiently reducing both u>o and u>h- 



In order to get more insight into the wave function we 
plot in figures 1-3 the (radial) pair correlation function 
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where the dimensionless coordinate £* is 

1/4 



r = 



V2l B 



(22) 



(23) 



This pair correlation function gives the probability of 
finding two particles at a given distance in the presence 
of a third particle. We find that, owing to the presence 
of the third particle, a shell structure appears with cer- 
tain maxima in the probability at special distances. The 
number of these maxima (shells) is defined by the quan- 
tum number n 2 and is equal to n 2 + 1. Therefore, the 
quantum number n 2 can be interpreted as the main char- 
acteristic of the number of resonant states or the number 
of nearest neighbors if one considers this as the onset of 
a liquid like behavior. 

This interpretation is supported by the following rea- 
soning. The pair correlation function starts at zero since 
the antisymmetrization due to Pauli-blocking induces the 
exchange correlation hole. A small £ expansion then 
yields [appendix A] 

g (\c\) = a\c\ 2m2 + (b+ cir i 2m2 )in 2 + o(in 3 ) 

_ r (a + 6)|r| 2 , m 2 = l 

- \ 6iri 2 + ain 2m2 , m 2 > 1 



+o(\en (24) 
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2. The pair correlation function <?(£*) versus £* for 
and different values of n 3 and m 3 . 



0.6 
0.4 




0.6 




0.6 
0.4 
0.2 


0.6 
0.4 



'\ n 2 =3 m 2 =1 n 3 =1 m 3 =0 




i i n 2 =3 m 2 =1 n 3 =1 m 3 =1 




- ' \ n 2 =3 m 2 =1 n 3 =3 m 3 =2 




n 2 =3 m 2 =3 n 3 =3 m 3 =3 





2.3 4 



2 . 3 



(iri) 



n 2 ! 



(n 2 + |m 2 |)! 



(4T'(r 2 ))' 



FIG. 3. The dependence of <?(£*) on space coordinate for 
112 = 3 and different values of ri3 and mz. 

The characteristic property of the Laughlin wave func- 
tion is that g(z) goes to zero as z 2m , whereas in the 
case of the Wigner crystal the pair-correlation func- 
tion goes to zero as 1 — exp (— az 2 ) cx z 2 . This was 
the reason for the superiority of the variational Laugh- 
lin wave function. In our direct computation, starting 
from the expression of antisymmetric wave function *S?Q , 
we obtain the Wigner crystal behavior. However, the 
Laughlin exchange correlation hole becomes transpar- 
ent if we use only the first diagonal term in | \I> | 2 cx 
C 2 |* ni , mi (r;w = 0)| 2 |*„ 2im2 (0| 2 |*„ 3 ,m3(??)| 2 which is 
the first term in (Al). We obtain the analytical result 



(25) 



Representing the Laguerre polynomials as L™(x) 



m 
k 



En 

ior of g dmg is obtained 



the following asymptotic behav- 



7 diag (|£*|) = 



m 2 



-r 2|m2 ' (i + o(e 2 )). (26) 



n 2 !m 2 ! 2 

We see that the typical z 2 !" 12 behavior of Laughlin states 
appears. Therefore, we consider the quantum number 
|m 2 | as the filling factor in analogy to Laughlin's ex- 
change correlation hole behavior. 

The approximation leading to this behavior has con- 
sisted in neglecting certain crossed terms. Here we should 
note that when considering the paircorrelation function 
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we design in principle two out of the three particles differ- 
ently than the third particle. In this case we don't need 
to antisymmctrize the wave function. One can therefore 
argue that all terms of Eq.(Al) except the first one do 
not count and the result (25) is exact. 

The function ,g dia s 1 in Eq. (25) is plotted in Figs. 1-3 
with dashed lines. We see that the places of the max- 
ima are well described by this approximation. Note here 
that the approximation g^s 1 indeed shows the Laughlin 
exchange correlation hole of £ 2 ™ 2 . 

In order to obtain analytical results for observables 
we can use all three diagonal elements |\[' < 3 dla s| 2 = 

C 2 ^! 2 (|* 2 | 2 |* 3 | 2 + |* 2 | 2 |* 3 | 2 + |<l 2 | 2 |^ 3 | 2 ). The 
advantage of this approximation is that one can com- 
pute the expectation values of observables analytically 
and get a fairly good approximation. 

We consider the mean distance of the particles 



(ri -r 2 f 



) = J d 2 Zd 2 rd 2 V e\^ QAias \ 2 
3C 2 l 2 

1 r (l + 2n 2 + |m 2 |) + (l + 2n 3 + |m 3 | 



(27) 



as well as the area spanned by the three particles, 
(i|(r 2 -n) x (r 3 -n)| 2 ) = J d 2 Zd 2 rd 2 r,^ x r,\ 2 \^ di ^\ 



\-[^k [(1 + 2n2 + |m2|)(1 + 2n3 + |m3 ' 



(28) 



We see that the mean distance and the mean area are 
determined by the quantum numbers 1 + 2n 2 + |m 2 | and 
1 + 2n 3 + |m 3 |. The mean distance is given by their sum 
while the mean area is given by their product. The non- 
diagonal terms give smaller and smaller contributions. 
In table I we compare the analytical approximate results 
(27) and (28) with the numerical value of the full wave 
function. The agreement especially for the mean distance 
suggests that this is quite a good approximation. 
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30 
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8.18 


30.50 
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27 


3332 


2.28 


8.64 


32.36 
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9.5 


33.75 


3133 


2.37 


8.15 


23.10 
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30 


3333 


3.70 


9.29 


41.88 




3 


10 


37.5 


TABLE I 


The normalization constant C, the dimension- 



less mean distance and the area for different quantum num- 
bers. The corresponding second lines give the analytical val- 
ues (27) and (28) respectively. 
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Although the states in high magnetic fields are char- 
acterized by spin alignment and are therefore described 
by the quartet states, we now give for completeness the 
other possible state, the doublet state. The wave func- 
tions corresponding to the doublet states can be con- 
structed according to Young's scheme [18]. For the state 
\S = 1/2, s z = 1/2) it is given in Fig. 4. The expressions 
corresponding to the configurations (a), (b), and (c) in 
Fig. 4 are obtained by antisymmetrizing the spin wave 
function with respect to the row indices and symmetriz- 
ing with respect to the column ones, 

1 1/2, l/2)„ = c (a 2 {a 1 (3 3 - a 3 /M + ai{a 2 /3 3 - a 3 f3 2 }) 

(29) 

and correspondingly for (b) and (c). It is possible to see 
that one particular spin wave function transforms into 
the other one when two particle indices are interchanged, 
forming a group. 
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FIG. 4. The particular Young diagrams contributing to the 
spin state \S = 1/2, s z — 1/2) are shown on the left with the 
corresponding contributions to the spatial wave function on 
the right. The indices in one row are understood to be sym- 
metrized and the indices in one column are antisymmetrized. 
The diagrams for \S = 1/2, s z — —1/2) are obtained by inter- 
changing a <-> [3. 

The spatial wave functions for each spin diagram are 
obtained by transposition of the Young tableaux of Fig. 4. 
Once again, the elements in a row must be symmetrized 
and the elements in a column antisymmetrized. The spa- 
tial wave function corresponding to (a) in Fig. 4, Eq. (29), 
reads 

tf £ (n, r 2 , r 3 ) = C^> numi (r; lu = 0) 

X {*n 2 ,m 2 {-Z)*n 3 ,mM + *n 2 ,m 2 (0*n 3 ,m 3 (v) 

(-|)*»s,rr,sW)}- ( 3 °) 

Denoting only the spatial indices, the coordinates £,77 
defined in Eq. (2) can be abbreviated as (1,23). The 

coordinate £,77 of (20) would correspond to (2,31), £,fj 
to (3,12), and —£,fj to (3,21) respectively. In this short 
notation the three orbital wave functions of figure (4) 
read 

*? = (2, 13) + (2, 31) -(1,23) -(1,32) 
*? = (1,32) + (1,23) -(3, 21) -(3, 12) 
? = (3, 21) + (3, 12) - (2, 31) - (2, 13). (31) 

Multiplying the spin and spatial components and 
adding all three particular expressions, we get the final 
wave function for the doublet state. The final result de- 
pends on the parity of m 2 . For m 2 odd it reads 

* D (r!,r 2 ,r3) -C7i* ril , mi (r;w =0) 

x{(/?ia 2 a3 + ai/3 2 a 3 - 2aia 2 /3 3 )*„ 2!m2 (£)*n 3 ,m 3 (v) 

+ (ai/3 2 a 3 + aia 2 /3 3 - 2/?ia2a 3 )*n 2 ,m 2 (0*n 3 ,m 3 (?7") 

+ {(3ia 2 a 3 + a 1 a 2 (3 3 - 2ai/3 2 a 3 )*„ 2 , m2 (f)*„ 3 , m3 (r))} (32) 

while for m 2 even the wave function is given by 

y D (r 1 ,r 2 ,r 3 ) = C7 2 *„ liTOl (r; uj = 0) 
x{(ai/3 2 a 3 - /3ia2a 3 )* n2 , m2 (0)*n3,m 3 (»?) 
+ (a 1 a 2 (i 3 - ai(3 2 a 3 )^ n2 . m2 (l)^ n3 . m3 (rf) 

+{(3ia 2 a 3 - aia2/33)*n2,m 2 (£)* n 3 ,m 3 

(33) 

It is easy to check that the expression for the proba- 
bility distribution \^ D \ 2 in the doublet state differs from 
the one for the quartet state by the coefficients in front 
of the crossed terms. Our estimate shows that the char- 
acter of quantization of the two-particle correlator, the 
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cluster area, mean distance for doublet states is qualita- 
tively similar to that for quartet states. Since the con- 
tributions of the crossed terms in the quartet states are 
smaller than those in the doublet states, the maxima in 
the pair correlation function are more pronounced in the 
former case. 



IV. CONCLUSIONS 

The investigation of the three-particle problem in the 
presence of magnetic fields is essential to understand 
the phenomenon of the fractional quantum Hall effect 
[2]. In this paper we have solved exactly the problem 
of three identical particles interacting via a spring po- 
tential in 2D space in the presence of an external mag- 
netic field, neglecting the Coulomb interactions. We have 
shown that the wave function acquires the Jastrow pref- 
actor even in the case of an attractive potential as a 
consequence of the anti-symmetric character of the to- 
tal wave function. Calculations of the pair correlation 
function, the area of the three-particle cluster and the 
mean distance show a quantization of these parameters. 
Particularly, the pair correlation function exhibits max- 
ima and minima controlled by a quantum number. 

We have also found that the exchange correlation hole 
at small distances shows the behavior known from the 
Wigner lattices rather than the Laughlin behavior if we 
start from the antisymmetrized wave function. However, 
neglecting certain terms based on the physical consider- 
ation that the antisymmetrization of the wave function 
fails when two particles are fixed for the pair correlation, 
we obtain exactly the Laughlin exchange correlation hole. 
The corresponding filling factor is identified with a quan- 
tum number appearing in the exact solution of the model 
solved here. 



APPENDIX A: SMALL DISTANCE EXPANSION 
OF THE PAIRCORRELATION FUNCTION 

The expansion of the paircorrelation function (22) for 
small £ can be performed straightforwardly with the help 
of MATHEMATICA. After trivial integrations about d 2 r 
we have to calculate 

g{\C\) = J dV(l* 2 | 2 |* 3 | 2 + 2|* 2 | 2 |* 3 | 2 

+2Re {2*2*3*2*3 + 2*2*3*2*3}) ( A1 ) 

where we have used symmetries in the variables and ab- 
breviated * 2 = *n 2 ,m 2 (0> *3 = *n 3 ,m 3 (^) and similar 
for the other combinations. Therefore we express the oc- 
curring variables (20) in the wave function in terms of 
the integration variables 



Iff 
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3M 2 + |£| 2 -2V3|77||£|cos(^ 



(A2) 



and similarly for fj as well as fj and £. Using the dimcn- 
sionless coordinates (23) in the wave function (14) and 
changing variables (A2), the integral (Al) is quite massy. 
However, expanding in small £* the angular dependence 
(f) v — <p£ drops out and we obtain the form (24) 

g(\€\) = a\C\ 2m2 + (b + c\C\ 2m2 )\C\ 2 + o(\C\ 3 ) (A3) 

where the remaining modulus integral d\r}\ can be per- 
formed for the constants a and c with the result 



a = 2 



(n 2 + m 2 ) 
m 2 ! 2 n 2 ! 



_ 2(n 2 + m 2 )! . , 
n 2 !m 2 ! 



n 2 



(1 + to 2 )! m 2 ! 



ZT7 • (A4) 



Therefore we see that the constant a is not vanishing 
which means that the Laughlin liquid behavior is not 
resumed if full antisymmetrization is pertained for the 
pair correlation function. 

For completeness let us also give the constant b which 
integral can be done only numerically [x = \r]*\ 2 ] 



b = 



^_ j^2m 3 gm2 — 1 



n 2 !n 3 ! 



4m 2 +m 3 {m 2 + n 2 y.{m 3 + n 3 )\' 

OO 

J dxx m2+m3 - 1 e- :c 


x|8(m 2 -3™ 3 ) 2 ^(|) 2 L™ 3 (i) 2 
St* t 

+12(m 2 -3m 3 )xL^( T )L™»( i ) 



,3x 
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